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Enhancement and Suppression of Four wave mixing in coupled semiconductor 

quantum dot-gold grating systems. 
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We have shown that Four wave mixing (FWM) processes of electromagnetic field modes of a 
grating can be controlled by the presence of interactions with a quantum dot or a molecule made 
by coupled quantum dots. By choosing the appropritae level spacing for the quantum emitter, one 
can either suppress or enhance the Four wave mixing process. We revel theoretically the underlying 
mechanism for this effect, (i) Suppression in FWM intensity occurs simply because induced Elec¬ 
tromagnetic Induced Transparency does not allow the excitation at converted FWM frequency, (ii) 
Enhancement emerges since FWM process can be brought to resonance. Path interference effect 
cancels the nonresonant frequency terms. Furthermore, we have also shown that in case of coupled 
quantum dots enhancement increases significantly as compared to the case of a single quantum dot. 


I. INTRODUCTION 

Quantum Plasmonics is an emerging area of research which involves the study of the optical properties of hybrid 
photonic structures incorporating both plasmonic nanostructures and quantum emitters [l], such as atoms, molecules 
and semiconductor quantum dots. These complex hybrid are active photonic structures and expected to enhance 
optical response significantly, for example modification of the linear susceptibility and the enhancement of 

nonlinear susceptibilities in several quantum systems with different level structures coupled to various plasmonic 
nanostrcutures j7Hll|. 

Four Wave Mixing is one of the above mentioned nonlinear process of light-matter interactions in which three 
incoming waves,indicated as uq, uq, W 3 in the material generate a fourth wave of frequency C 04 pjj]. Assuming that 
the three incident waves have frequencies in the visible or near-infrared range, the incoming electric fields E (uq) 
(with i = 1 , 2 ,3) interact with the material’s electrons to induce a nonlinear polarization P^ (uq) in the illuminated 
volume. The magnitude of the polarization is determined by the strength of the incident fields and the efficiency 
with which the material can be polarized [l3|. The latter is indicated with the third-order nonlinear susceptibility 
a measure of the material’s response to the incoming fields. Four wave mixing (FWM) has found numerous 
practical applications, including: optical processing; nonlinear imaging; real-time holography and phase-conjugate 
optics; phase-sensitive amplification; and entangled photon pair production [l3 |. 

In several recent studies, the modification of x ® (FWM process) susceptibility in a quantum dot system coupled 
to spherical nanoparticle has been investigated when the hybrid structures interacts with a weak probe field and 
a strong pump field All these works have shown for different distance between the quantum dot and 

the metal nanoparticle the susceptibility can be either enhanced or strongly suppressed. In addition, bistable 
behavior has been also reported in these kind of systems EMI- 

Here, we propose a method for increasing the efficiency of FWM processes by exploiting gold grating 51111 - Narrow 
peaks are observed in the transmission spectra of p-polarized light passing through a thin gold film that is coated 
on the surface of a transparent diffraction grating. The spectral position and intensity of these peaks can be tuned 
over a wide range of wavelengths by simple rotation of the grating p5j ). The wavelengths where these transmission 
peaks are observed correspond to conditions where surface plasmon resonance occurs at the gold-air interface. Light 
diffracted by the grating couples with surface plasmons in the metal film to satisfy the resonant condition, resulting 
in enhanced light transmission through the film. 

The paper is organized as follows. In Section II, we describe the FWM Process in the coupled system of gold grating 
with a quantum oscillator. In the same section, we introduce the Hamiltonian for hybrid system. FWM process is also 
included in the second quantized Hamiltonian. We derive the equations of motion for the system using the density 
matrix formalism for the quantized quantum oscillator. We use phenomenological way to include damping of gold 
grating modes as well as quantum emitter. In Section II B, We demonstrated that FWM process can be suppressed 
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for (w eff =2w — a/). In Section II C, we present a contrary effect where FMW process can be enhanced. This is due 
to cancellation of non resonant terms in denominator. In Section III, we further investigate the case of gold grating 
coupled to two quantum emitters (both quantum emitters coupled to each other also) simultaneously. We conclude 
our results in Section IV. 


II. HAMILTONIAN FOR FOUR WAVE MIXING 

The total Hamiltonian H for the described system can be written as Sum of the energy of the Quantum Oscillator 
Hq,( In our case we have taken a QD of energy levels |e) and |c/)) , enegy of the elctromagnetic modes of Gold Grating 
(ai, a 2 , a 3 ) for a particular angle of incidence of Pump lasers Hgrating, the interaction of the Quantum Oscillator with 
the Grating Modes H int 


H 0 = hu e \e) (e| + hu g \g) (g| (1) 

H grating = hbj\a\a\ + h,L0 2 d\d 2 + tuj 3 dla 3 (2) 

Hint = h (/ 4 Ip) (e| + / *a 3 \e) (g|) (3) 

Here, we have considered that level spacing of the QD is only resonant to a 3 mode (i.e. u> eg ~ uj 3 ). as well as the 
energy transferred by the pump source ui and a/. 

H P = ih (4 c P e~ iut - cl 1 e* P e iujt ^j + ih (o£ - a 2 eJJe 4 "'*) (4) 

H FW m = S.X (2) ( a\a\a\ + a\ 2 a 2 a 3 ^ (5) 

For the Process of (uj 3 = 2uq — u> 2 ) as mentioned in PRL 103, 266802. 


In Eq. (1), hco e (hu! g ) is the excited (ground) state energy of the Quantum Oscillator. States \e) and \g) corresponds 
to excited and ground levels of the Quantum Oscillator respectively. ( 01 , 02 , 03 ) are the Gold Grating modes at a 
particulat angle of incidence 0 = 5°. f is the coupling matrix element between the field of grating mode and the 
Quantum Oscillator. Eq. (4) describes the driving the Electromagnetic field modes of grating (ai and a 2 ) with e~ lut 
and e~ lul t respectively. 

Eq.(5) describes where the Four wave mixing takes place in which ai mode contributes two photons and a 2 mode 
single photons in the process. 


A. Heisenberg Equations of Motion 


We use the commutation relations 


ihSo =\o,h 
dt L J 

( 6 ) 

for deriving equations of motions. After obtaining the dynamics in the quantum approach, 
classical expectation values (011,0:2,0:3) ■ We also introduce the decay rates for (01,02,03). 
treated within the density matrix approach. The equations of motion take the form 

we carry (di,d 2 ,a 3 ) to 
. Quantum Oscillator is 

di = (—iu) 1 — 71) 01 — 2ix^ala 2 a 3 + e p e~ lult 

(7a) 

02 = (—ioj 2 — 72) 02 — i X (3) o*o? + c p e~ lu 4 

(7b) 

d 3 = (~ itjJ 3 - 73) o 3 - *X (3 ) 02Qi - ifpge 

(7c) 
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Pge — ( i^eg 7ep) Pge 7 ifo 3 (Pee Pgg) 


(7d) 


Pee — ^ieePee 7 if (oq/3 g e 03Peg) 


(7e) 


where 71,72,73 are the damping rates of the electromagnetic modes of the gold grating ( 01 , 02 , 03 ) ■ 7 ee and 7 eg = 
7 ee /2 are the diagonal and off-diagonal elements of the quantum oscillator respectively. The constraints of the 
conservation probability p ee 7 p gg = 1 accompanies above set of equations. 

Besides the time-evolution simulations, one may gain the understanding by seeking solutions of the following form. 
For long time behavior we take solutions of the form 

a\(t) = a\e~' lult , a 2 it) = 02 e~ tul ' t , a^ft) = (Condition of Four wave mixing process), p ge = 

Pgee-^ l ^ 2u ~ w here we have considered that level spacing of the QD is only resonant to 0,3 mode (i.e. ui eg ~ W3), 
Pee — Pee- 

Inserting the solutions in above set of equations (7a-7e) we have for long time behavior 


[i (uq — uj) 7 71] di + ot\ 6 t 2 (X 3 = < 

[i (u 2 - u/) + 72] a,2 7 *x (3) d3di = e(, 

[i (W 3 7 <*/ - 2 uj) + 73 ] d 3 + ix^a^al = -ifp, 

[i (u) eg 7 J - 2u) + 7eg] Pge = ifo 3 ( Pee - p gg ) 
7 eePee = if [,^ 3 Pge O^Peg) 

Using equations ( 8 c) and ( 8 d), we obtain the steady state value for 0:3 as follows 


ge 


03 — 


*X (3) d 2 «i 


:i( u .,+ 5 LL)+ 7 . e ] _ I* ( w 3 + w' - 2 w) + 73] 


( 8 a) 

( 8 b) 

( 8 c) 

( 8 d) 

( 8 e) 

( 9 ) 


Where y = ( p ee — p gg )is the steady state value of the population inversion. If the quantum oscillator is tuned around 
u) eg = 2 u> — u/, oq can be suppressed. 


B. Suppression of the Four wave Mixing Process 

We can see from Eq. (9) that \.f \ 2 y/jeg can attain huge values on resonance uj eg =2 uj — u/as well as linewidth 
of the quantum oscillator 7 eg is very small compared to the all other frequencies. If / ^ 0, the largeness of the 
l /| 2 y/leg term dominates the denominator. This results in the suppression of the generation of the FWM mode S 3 
in our model Hamiltonian system. In Fig.l we have shown that FWM process can be suppressed very effectively by 
coupling to gold grating to quantum oscillator. We have time evolve Eqs. (7 a — 7e) to obtain steady state values for 
the FWM intensity. 

i ( 3 ) q :* q : 2 

Without the presence of quantum oscillator, the FWM would be maximum S 3 = w h en th e FWM mode is on 

resonance (u !3 = 2uj — u>'). In Fig. 1 we observe that even at the presence of this resonance condition (w 3 = 2 uj — uj'), 
EIT suppresses the FWM by 10 order of magnitude. Furthermore, in this case population of excited level of quantum 
oscillator is maximum at this point as shown in Fig.2 as well as population inversion is approximately y = ( p ee — p gg ) « 

1 


C. Enhancement of Four Wave Mixing Process 


Similar to suppression phenomena, the interference effects can be arranged in such a way that FWM process can 

| r I 2 

be carried closer to resonance. In the denominator of Eq.(9), the imaginary part of the first term — -hJ'- 2 u )+7 —f 
can be arranged to cancel the i (uj 3 + uj' — 2uj) factor in the second term of the denominator. This gives the condition 


l/| 2 V ( w eg + Ul' — 2Uj) + (u >3 + U)' 


2uj) (uj eg 7 uj — 2tu) 7 7 , 


= 0 


( 10 ) 
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FIG. 1: Suppression of the FWM intensity to the 03 gold grating mode from the a 2 and a\ mode. Even at the 
presence of the resonant FWM condition, w 1 = w = 1.0, 002 = 00 ' = 0.5 and (w >3 = 2w — u/), the presence of quantum 
oscillator prevents to take place of the FWM process. EIT does not allow the FWM process. The resonant FWM 
conversion is represented by unity in figure. When (ui eg = 2w — ui'), the FWM intensity even can be suppressed by 
10 orders of magnitude with respect to resonant value. Decay rates for our numerical simulations are 
7 i = 72 = 73 = O.Olw and 7 eff = O.OOOOlw. We have taken %^ = O.OOOOlw and / = O.lw 



FIG. 2: Enhancement of population in excited level of the quantum oscillator coupled to gold grating at the 
resonance condition ( oj eg = 2w — u/). We can see the population in excited level is maximum at this condition unlike 
FWM Intensity | S 3 | 2 shown in Fig. 1. All other parameters for numerical simulation remain same like in Fig 1. 


Eq.(10) has two roots. 


(veg 2) +oj' - 2w) = 


\f\ 2 y 


T 


l/IV 


( 0 J 3 + ui' — 2 w) y -\- ui' — 2 cjj)‘ 


-4 7 | 


( 11 ) 


The hrst smaller root u>i g ~ 2 u — uj' is not very useful for FWM enhancement, as it enhance the real part of the term 
— + jJ}_ 2 U ) +/y —j to rapidly diverge as we have seen in suppression condition for FWM, whereas wig* minimizes the 
absolute value of the denominator of Eq. (9) that gives enhancement of FWM process. For the case of suppression of 
FWM, one can safely use the approximation y « — 1 because excitations are suppressed in the hybrid system p ee « 0 
and this leads to y =(p ee — p gg ) ~ —1. However, in case of FWM enhancement, one can not approximate y s=s —1. 

Nevertheless, Eq.(ll) still serves at least a guess value for the order of wig\ where FWM enhancement arises. 
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FIG. 3: The enhancement of the FWM process. The FWM mode <23 is far-off resonant to the FWM condition 
(023 = 1.85w). The FWM process can be carried closer to resonance by arranging the quantum level spacing to 
( uj eg ss 1.52w). The conversion is enhanced nearly 80 times compared to off-resonant process. The conversion for off 
resonant process (/ = 0) is represented by unity in figure. For {u> eg = 2w — a/), FWM process is suppressed similar 
to Fig.l. Decay rates of the grating modes are taken as 71 = 72 = 73 = O.Olw. We use x^ = O.OOOOlw and / = O.lw 

for numerical simulations. 


III. TWO COUPLED QUANTUM DOTS 


In case of two coupled QDs we have total Hamiltonian as Follows. 


Ho = Mo l e i) <ei| + Ku ( eq l e 2 > (e 2 | 


Hgrating ~ + hid 2^ 2 ^2 + hu 3 d 3 d 3 


Hint = h Ifl'i) (ei| + /i *03 \ei) (fli|) + (f 2 a\ |p 2 > (e 2 | + / 2 *a 3 |e 2 ) ( 52 1) (14) 

Here, we have considered that level spacing of both the QDs is only resonant to 0,3 mode (i.e. u>lg* ~ W 3 , Weg'* ~ W 3 ). 

Hqe—qe =h\g (|e 2 ) (g 2 | ® Ifl'i) (ei|) + fl* (|ei) (fli| ® |fl 2 ) (e 2 |)] (15) 

as well as energy transferred by Pump Source with e~ lult and e~ lu t respectively. 

Hp = ih e P e~ iut - ax e* P e iujt ^j + ih (aj e^e "^' 4 - a 2 (16) 


Hfwm = ?b\' (3) ( a\a\a\ + a\ 2 a 2 a 3 ) 


is the Four Wave Mixing Hamiltonian for the Process of (oj 3 = 2uq — w 2 ) as mentioned in PRL 103, 266802. 

By using the commutation relation Eq.( 6 ) as well as proceeding like the same way for a single QD case, we get 
following equations for the case of coupled QDs. 

ax = (— iuj\ — 71 ) a.\ — a*a 2 a 3 + e p e~ lult (18a) 


Oil = {-iu 2 - 72 ) ol 2 - *x (3) a 3 aj + e' e * 


d 3 = (-iw 3 - 73) «3 - ix (3) a 2 o;? - ifip$ - if2P { g J 


( 18 c) 
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nW = 


( 1 ) 


Pge = [ ~^ea ~ 7, 


eg 


v(1) ) pfe + i/l*«3 (/4e } - Pgg) + W* {p ( ee ~ Pgg } ) 


'e</ / rge 


p (1) 1p (2) 


pie = - 7^) P { g} + */ 2 *«3 (ple } - P$) + ig (pl 2) - P$) p! 


( 2 )_ rt ( 2 )U(l) 


99 / 


pll } = -7llVll ) + i (jia%p { g} - /ra 3 p^) + i (gpfJpQ ~ 9*P { eg P$ 

P^ee = ~1ee Pee + < (/^ ~ + * (flVM ~ *»«) 


(18d) 

(18e) 

(18f) 

(18g) 


where 71 , 72,73 are the damping rates of the electromagnetic modes of the gold grating ( 01 , 02 , 03 ). 7 i)\ 'yif and 
leg ' 1 = 7 ee^/ 2 , 7 ^ = 7 ee^/ 2 , are the diagonal and off-diagonal decay rates of the first and second quantum emitter 
respectively. The constraints of the conservation probability p£) + p'gg = 1 and p'e) + pig 1 = 1 accompanies above 
set of Eqs.(18a-18g). 

In our simulation for enhancement process of FWM, we time evolve Eqs. (18a-18g) numerically to obtain the long 
time behaviors of pi)) , pi))) , pi)) , pie'*, 01 , o 2 and 03 . We determine the values to where they converge when the drive is 
on for long enough times. We perform this simulations for different parameter sets ^/i, / 2 , g, u4g\ a4g\ 7 eg\ 7 eg^ with 
the initial condition pi)) (t = 0) = p)) (t = 0) = 0, PgV (t = 0) = Pg? (t = 0) = 0, a\ (0) = 0, a 2 (0) = 0, a .3 (0) = 0. 

Besides the time-evolution simulations, one may gain the understanding by seeking the solutions of the following 
form: 

a\{t) = aie~ lu>t , a 2 (t) = d 2 e - “‘ J,t , a 3 (t) = a 3 e _I ( 2 “ _ “ )* (Condition of Four wave mixing process), pi)) = 

~(1) (2) ~(2) —i( 2 uj—uj')t 

Pge £ V ' , pge = Pge € V ) . 

Here we have considered that level spacing of both the QDs is only resonant to S 3 mode (i.e. lu))) ~ 103 , 10 )) ~ w 3 ), 
pl^ (i) = pi)) and pi 2) (t) = pi). 

Inserting the solutions in above set of Eqs. (18a-18g) we have the following closed set of equations for the steady state 
dynamics 


[i (wi — u>) + 71 ] di + 2 ix^a\a 2 a 3 = e p (19a) 

[i (w 2 — to') + 72 ] d 2 + ix^dgQq = e' p (19b) 

[* (w 3 + u} 1 - 2w) + 73 ] 03 + *X ( 3 ) d 2 b l = -ifip)) -*/ 2 Pg 2) (19c) 

i (4j + w' - 2w) + 7 gg } p),) = ifi& 3 y 1 + ig*yipf) (19d) 

i (u))) + J - 2w) + 7 ^ 4? = */ 2 d32/2 + igyipf} (19e) 

7ie ) pii ) = i (/ldjiPge _ f* & 3P ( ej) + i (sPigPge “ 9* P { )) 0f)) (19f) 

7i 2 ) pi 2) = * (fiOllpf} - f^&spi))) + i (p*pig ) Pg 2) - PPigPge) ( 19 §) 


x(2) r 


where aq, d 2 , 0 I 3 , pge*, pg 2 \ pie'* and pip*are constants independent of time, yi = ^pie — PgPj are the population 
inversion (* = 1,2) for both QDs. 

Using Eqs.(19d) and (19e) in Eq.(19c), we obtain the steady state value for 0:3 as follows. 


013 — 


*x (3) (Pip 2 + yiy 2 |p| 2 ) 


(yi I/ 1 I 2 fc + 92 1/ 2 | 2 0i ) + iym (/i/ 2 *p* + /1/25) - £3 (/?i/ 3 2 + ym |p| 2 ) 




( 20 ) 


where the short hand notations are £\ = [i (wi — ui) + 71] , £2 = [i (w 2 — w') + 72] , £3 = [i (0J3 + u>' — 2 w) + 73] 

(2) _ 


and 0i = i (u>i)) + to' - 2 and 02 = i (uiif + u' - 2 w^ + 7 ^ 
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A. Super enhancement of FWM process 

1. Single QD case 

In case of a single QD coupled to the gold grating, f 2 =g = g*= 0 and f± = /, we get the steady state value of 
d 3 from Eq.(20)as 


0:3 


iX^alal 


i(u) eg +ll'-2uj)+~/eg] ~ ^ ^ + ~ ^ + 


( 21 ) 


I ci2 

which coincides exactly with Eq.(9) where the imaginary part of the first term —j can be arranged to 

cancel the i (w 3 + a/ — 2 oj) factor in the second term of the denominator and this gives enhancement of FWM as also 
discussed in previous section also. 


2. Coupled QDs case 

As compared to single QD case, the denominator of Eq.(20) can (in principal) be arranged down to very low 
values in order to enhance d 3 to much higher values. In this case, denominator has 3 complex (/1 , f' 2 , g) and 2 real 
,bJejj j parameters which can be tuned independently. 

We obtain nearly 1200 times enhancement by comparing the steady state values of |ck 3 | 2 that is the intensity of 
FWM process calculated from time evolution of Eqs.(18a — 18g) for the chosen set of parameters as shown in Fig.(4). 
Here, frequency of second QD is kept constant and first one is varying. For the decay rates of grating modes in between 
0.01, we get enhancement in FWM Intensity around 1200-1600 times as compared to case of single QD discussed in 
previous section. 


1600 

1400 

C\J_ 

JO 1200 
o 

1000 

>* 

c 800 

a> 

600 

o 

x 400 
200 


FIG. 4: The enhancement of the FWM process in case of coupled QDs. The FWM mode a 3 is far-off resonant to the 
FWM condition (w 3 = 1.90a;). The FWM process can be carried closer to resonance by arranging the quantum level 

spacing of first QD to ~ 1.5732u;^, while second QD (uJef = 1.5810u;^ being fixed . The conversion is 
enhanced nearly by 1600 times. Decay rates of the grating modes are taken as 71 = 72 = 73 = 0.01a;. We use 
X (2) = 0.00001a; and h = f 2 = 0.1909w, g = (0.1000 + 0.0101*) 077 ^ = 7 ^ = 0.00001a;, 7 $ = 7 ^/2, 7 ^ = 7^/2 

for our numerical simulations. 
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IV. CONCLUSION 


It is well demonstrated that the presence of a quantum emitter with a smaller decay rate changes the optical response 
of coupled grating dramatically. Due to the destructive interference of the (hybridized) absorption paths, Four wave 
mixing(FWM) process can be suppressed at the resonance frequency of the quantum emitter. We demonstrate that a 
similar path interference effect can be adopted to both suppress and enhance the nonlinear Four wave mixing processes 
(FWM) in a grating surface. A quantum emitter is coupled with the electromagnetic modes of a gold grating. We 
found that the FWM process can be suppressed over 10 orders of magnitude. Such an suppression can be achieved by 
carefully choosing the coupling strengths and the energy level spacing for quantum emitters. When (w es = 2w — ui'), 
the FWM intensity can be suppressed by several order of magnitude with respect to resonant value. On the other 
hand, the similar interference effects can be also used to enhance the nonlinear FWM intensity. The level spacing 
of the single quantum emitter can be arranged so that the nonresonant terms get canceled. In case of two coupled 
quantum emitters by arranging energy level spacing for quantum emitters in the same way like single quantum emitter, 
we have enhancement in FWM intensity upto the order of 10 3 . 
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